DFTUZ/96/3 



TESTING LOGARITHMIC VIOLATIONS TO SCALING 
IN STRONGLY COUPLED QED 



V. Azcoiti", G. Di Carlo^ A. Galante^'^ A.F. Grillo"', V. Laliena" and C.E. Piedrafita'' 



" Departamento de Fisica Tedrica, Facultad de Ciencias, Universidad de Zaragoza, 

50009 Zaragoza (Spain). 

^ Istituto Nazionale di Fisica Nucleare, Laboratori Nazionali di Frascati, 
P.O.B. 13 - Frascati OOOU (Italy). 

^ Dipartimento di Fisica dell'Universitd delVAquila, L'Aquila 67100 (Italy) 

Istituto Nazionale di Fisica Nucleare, Laboratori Nazionali del Gran Sasso, 
Assergi (L'Aqmla) 67010 (Italy). 



ABSTRACT 

Using very precise measurements of the critical couplings for the chiral 
transition of non compact QED4 with up to 8 flavours, we analyse the be- 
haviour of the order parameter at the critical point using the equation of state 
of a logarithmically improved scalar mean field theory, that of the Nambu- 
Jona Lasinio theory and a pure power law. The first case is definitively 
excluded by the numerical data. The stability of the fits for the last two 
cases, as well as the behaviour with the number of flavours of the exponent 
of the logarithmic violations to the scaling favour clearly a pure power law 
scaling with non mean field exponents. 
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Non perturbative renormalizable gauge theories are among the most stim- 
ulating and less known subjects in theoretical physics. It is not known if it 
is possible to construct such a theory in four dimensions and the standard 
prejudice is that only asymptotically free gauge theories with gaussian fixed 
points are "good" field theories. 

Non compact QED^^ or its generalized version, the gauged Nambu-Jona 
Lasinio {GNJL) model, are good candidates to analyze this problem. Both 
models have the common feature that they couple fermions strongly enough 
to produce fermion condensates, and therefore a phase, where chiral symme- 
try is spontaneously broken, appears at sufficiently strong gauge coupling. 
Since composite scalars are present in the spectrum of these models, the ex- 
istence of a non trivial continuum limit is strongly related, as discussed in 
[H], with the balance between the fermion attraction due to the interaction 
and the zero point repulsion due to the kinetic energy. If the short distance 
attraction is too strong, composite scalars with vanishing physical size ap- 
pear in the spectrum, thus giving rise to a non interacting field theory. Even 
if no rigorous proof exists, this seems to be the case in the Nambu-Jona 
Lasinio (NJL) model where mean field exponents with logarithmic viola- 
tions to scaling driving to a vanishing renormalized coupling are expected 

i- 

In a recent investigation of the gauged NJL model we have found 
strong evidence supporting the fact that the gauge interaction can change 
qualitatively the trivial scenario of the NJL model. Stimulated by these 
results and following our investigation on non compact QED4 started several 
years ago 0], we want to report in this letter some new results concerning 
the critical behaviour of this last model. 

Our analysis is based on very precise determinations of the critical cou- 
plings obtained from the computation of the susceptibilities in the Coulomb 
phase on 10^, 12^ and 14^ lattices. From these results we analyse the be- 
haviour of the order parameter at the critical point and for several number 
of fiavours by fitting it with three different equations of state (EOS): i) The 
EOS of a logarithmically improved scalar mean field theory, ii) a power law 
scaling without logarithmic violations and iii) a NJL model-like EOS. 

The case i), as will be shown in what follows, is definitively excluded by 
the numerical data. This is not surprising at all since as suggested in [|1| 
and corroborated in 0, triviality in a theory with composite scalars should 
manifest in a different way than in a theory with fundamental scalars. In 
fact the logarithmic violations to scaling in the first case, as in the NJL 
model, are expected to have an effective 6 exponent less than 3 whereas in 
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the second case an effective 6 larger than 3 is obtained. 

The difficult task is to distinguish from the numerical data between ii) 
(power law scaling) and iii) (four fermion EOS) since in both cases the ef- 
fective 6 is less than three. Let us anticipate that even if the behavior of 
the order parameter at the critical point is well fitted by both equations 
of state, the stability of the fits with lattice size and the behavior of the 
exponent which controls the logarithmic violations to scaling with the num- 
ber of flavours strongly favour a pure power law scaling against logarithmic 
violations to mean field. 

The main ingredient of this analysis is the very precise determination of 
the critical coupling f3c in QED^^ which follows from the computation of the 
susceptibility in the Coulomb phase of this model. A detailed explanation of 
this computation for lattices up to 10'^ can be found in p. 

Let us recall the two essential ingredients we have used in this analysis. 

First, exploiting the potentialities of the Microcanonical Fermionic Aver- 
age [MFA) approach to simulate dynamical fermions in gauge theories 
we have done simulations in the chiral limit at exactly zero fermion mass, 
which allows to overcome the ambiguities in the fermion mass extrapolations. 

Second, we have also exploited the fact that in the Coulomb phase of this 
model, characterized by a non degenerate ground state, we can interchange 
the chiral limit with the thermodynamical limit. Notice that the reason why 
this exchange is not allowed in general cases {QCD simulations, broken phase 
of QED, etc.) is the degeneration of the vacuum state. If we are in a broken 
phase where many degenerate vacua connected by symmetry transformations 
coexist, a permutation of the chiral and thermodynamic limits would imply a 
path integral over all the Gibbs state. Then several unpleasant facts like the 
violation of the cluster property for correlation functions, vanishing values 
for the order parameter or wrong values for the computed susceptibilities 
will appear. However if calculations are done in the symmetric phase char- 
acterized by a non degenerate ground state, there are not, in general cases, 
physical reasons to prevent from doing such a permutation. 

The most general expression for the vacuum expectation value of any 
operator O, after integration of the Grassmann variables, can be written as 
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where 
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n{E)= j[dA,]5{WE-SG[A^]) 



(2) 



is the density of states at fixed energy E and S^jj-{E, m) in (1) is the fermion 
effective action defined as 



O det A means the mean value of the product of the operator O times the 
fermionic determinant, computed over gauge field configurations at fixed en- 
ergy E. 

Since we are interested here in the computation of susceptibilities in the 
chiral limit, we will use the previous expression for the particular cases in 
which O is the longitudinal or transverse susceptibility. 

The longitudinal and transverse susceptibilities in the Coulomb phase are 
equal except a sign. They can be computed by taking for the operator O the 
expression 



where the sum in (4) runs over all positive eigenvalues of the massless Dirac 
operator. 

The transverse susceptibility xt in the broken phase diverges always be- 
cause of the Goldstone boson. The longitudinal susceptibility xl on the other 
hand, can not be computed in the broken phase using equation (4) since in 
this phase the ground state is not invariant under chiral transformations. 

We refer the reader interested to details on the computation of suscepti- 
bilities to 0. Let us recall that using the MFA approach, computations at 
several number of flavours Nf can be done without extra computer cost [^]. 
In table I we report the critical couplings extracted from the susceptibilities 
in lO'^, 12^ and 14^ lattices at Nf = 1,2,3,4,8. The most striking fact of 
this table is the small statistical errors of the critical couplings which follows 
from the high quality of the susceptibility fits, as reported in 

We have fitted the behaviour of the chiral condensate ( XX ) against the 
bare fermion mass m at the critical values of jS reported in Table I using 
different equations of state. In Fig. 1 we plot (xx)^ /^og{xx) against m 
in the mass interval (0.005,0.1) for the four flavour theory and in a 14^ 
lattice. Were the critical behaviour of this model described by a gaussian 
fixed point as in a logarithmically improved scalar mean field theory, the 
points in the small mass region of this figure should be well fitted by a 



log detA{E, m) 



(3) 




(4) 
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straight line crossing the origin. The sohd hne in this figure is the best fit 
obtained under the previous assumption. The very bad quahty of this fit 
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(doT ~ 1267, m < 0.1, = 358, m < 0.05) disproves definitively such 
a possibility. Similar results have been obtained in smaller lattices and for 
different fiavour numbers. 

Having shown that the critical behaviour of non compact QED can not be 
described by a logarithmically improved scalar mean field theory, we will try 
the same kind of fit for the other two reasonable possibilities, i.e. pure power 
law behaviour and mean field with logarithmic violations a la Nambu-Jona 
Lasinio. 

Fig. 2 is a plot of the same data reported in Fig. 1 but in the ordinate axis 
we plot ( XX )^'^^- The solid line is the best fit of all the points at small m with 
a straight line crossing the origin. The high quality of this fit is corroborated 
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by the value = 0.38, a value which is stable until masses of the order 
of 0.08 and increases slowly when masses larger than 0.08 are included in 
the fit. The value 2.80 of the 6 exponent has been chosen as the best one 
for the linear fit (see Table II and the discussion on the determination of 6 
at the end). From the results reported in Fig. 2 we conclude that a pure 
power with 6 = 2.80 describes with high accuracy the behaviour of the chiral 
condensate at the critical coupling. 

Our last plot for the chiral condensate at the critical coupling is reported 
in Fig. 3. For this plot we use a equation of state a la Nambu-Jona Lasinio 

(xx)W(7^)=Cm (5) 

where C in (5) is a constant and the exponent p of the logarithmic violations 
to scaling is left as a parameter of the fit. Recall that p = 1 in the large Nf 
limit and that this result does not changes after taking into account the 

M\ and B corrections, this suggesting that a value of p different from 1 is 

not very reliable. Notwithstanding that we decided to left p as a parameter 
of the fit for two reasons: i) our results for the chiral condensate does not 
support a fit with equation (1) and p = 1 and ii) there is no rigorous proof 
that p does not changes with Nf. 

The best fit of our results in the four fiavour model with equation (5) is 
obtained for p = 0.28. In Fig. 3 we plot {xx)^^og^''^^{j^^) against m. The 
solid line in this figure is a fit of all the points at small m with a straight line 
crossing the origin. Again now, as in the case of the power law behaviour, 

2 

we get a high quality fit of these points (^^^ = 0.89). 
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At first sight it seems difficult to decide between the last two cases (power 
law and Nambu-Jona Lasinio). However we can get some insight on the 
reliability of the last two fits by analyzing their stability with the lattice size 
as well as the fiavour dependence of the p exponent in equation (5). 

Table II contains a summary of the results we obtain for the 6 and p 
exponents for the three lattice sizes 10,12,14 and Nf = 1,2,3,4,6,8. The 
values of 6 reported in this table are obtained fitting the chiral condensate 
with a pure power law (case ii previously, corresponding to non mean-field 
exponents). Using instead a Nambu-Jona Lasinio like EOS, where 6 = 3, 
we obtain values of the p exponent of equation (5), which we also report in 
Table II. The errors in this Table take into account both, the errors of the 
fits and the error in the determination of the critical couplings (see Table 
I). The high precision of the exponents reported in Table II follows from the 
very small errors in the critical couplings of Table I. 

Looking at the results reported in Table II we can discuss both, the 
stability of the results with the lattice size and the fiavour dependence of p. 
First notice that the values of 6 for the three lattice sizes at each value of 
Nf are always compatible whereas this does not hold for the values of p. We 
conclude that the pure power law fits are much more stable with lattice size 
than A^JL-like fits. 

The second important fact that can be observed in this Table is that the 
value of the p exponent, contrary to expectations based on the expansion 
of the NJL model, not only is different from 1 but decreases when Nf in- 
creases and therefore goes away from the expected p = 1 at large Nf. As well 



known 0,|T^ the chiral transition of noncompact QED changes from second 



to ffrst order at large Nf and therefore we can not extrapolate our results for 
p in Table II to Nf = oo. Notwithstanding that, Kim, Kocic and Kogut have 
shown in ^ how numerical simulations of the NJL model with discrete Z2 
symmetry at Nf = 12 reproduce very well the logarithmic violations to mean 
ffeld scaling given by the EOS of the model at Nf = 00. Therefore, were 
the critical behaviour of non compact QED described by the NJL model, 
we would expect a value of p ~ 1, at least at Nf = 8. 

In conclusion, we believe that the features previously discussed strongly 
favour a pure power law scaling with non mean ffeld exponents against log- 
arithmic violations to mean ffeld a la Nambu-Jona Lasinio. This is a very 
interesting and important improvement with respect to previous work on this 
subject and it has been possible as a consequence of the very precise 



measurements of the critical couplings reported in Table I. 
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Figure captions 

Figure 1. Numerical results for {xx)^/^og{xx) against m for the four 
flavour theory and in a 14^ lattice. The solid line in this figure is the best 
linear fit crossing the origin. 

Figure 2. Plot of ( XX )^ '^° against m. The solid line is the best fit of all 
the points at small m with a straight line crossing the origin. 

Figure 3. ( XX against m. The solid line is a fit of all the 

points at small m with a straight line crossing the origin. 
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Table captions 



Table I Critical couplings extracted from the susceptibilities in 10^, 12^ 
and 14'' lattices at Nf = 1, 2, 3, 4, 8. 

Table II Results for the 6 and p exponents for the three lattice sizes 
10,12,14 and Nf = 1,2,3,4,6,8. The values of 5 reported in this table 
correspond to the results for the fits of the chiral condensate with a pure 
power law whereas in the case of the Nambu-Jona Lasinio like EOS, where 
S — 3, we report the value of the p exponent in equation (5). 
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